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ABSTRACT 

We present the first multi-color view of the scattered light disk of the Herbig Ae star HD 163296, 
based on coronagraphic observations from the Hubble Space Telescope Advanced Camera for Surveys 
(ACS). Radial profile fits of the surface brightness along the disk's semi-major axis indicates that the 
disk is not continuously flared, and extends to ~540 AU. The disk's color (V-I)=l.l at a radial distance 
of 3'.'5 is redder than the observed stellar color (V-I)=0.15. This red disk color might be indicative of 
cither an evolution in the grain size distribution (i.e. grain growth) and/or composition, both of which 
would be consistent with the observed non-flared geometry of the outer disk. We also identify a single 
ansa morphological structure in our F435W ACS data, which is absent from earlier epoch F606W and 
F814W ACS data, but corresponds to one of the two ansa observed in archival HST STIS coronagraphic 
data. 

Following transformation to similar band-passes, we find that the scattered light disk of HD 163296 
is 1 mag arcsec -2 fainter at 3'.'5 in the STIS data than in the ACS data. Moreover, variations are 
seen in (i) the visibility of the ansa(e) structures, in (ii) the relative surface brightness of the ansa(e) 
structures, and in (iii) the (known) intrinsic polarization of the system. These results indicate that the 
scattered light from the HD 163296 disk is variable. We speculate that the inner disk wall, which Sitko 
et al. suggests has a variable scale height as diagnosed by near-IR SED variability, induces variable 
self-shadowing of the outer disk. We further speculate that the observed surface brightness variability 
of the ansa(e) structures may indicate that the inner disk wall is azimuthally asymmetric. 



Subject headings: circumstellar matter — stars: individual (HD 163296) 
formation — planetary systems: protoplanetary disks 



planetary systems: 



1. INTRODUCTION 

Investigations of the nature and evolution of dust grains 
in proto-planetary systems are motivated in part by our 
desire to understand the birth and evolution of planetary 
bodies, which originate from these systems. It is well es- 
tablished that Herbig Ae/Be stars ( Herbig|1960 1 are inter- 
mediate mass pre-main-sequence stars, analogous to the 
more familiar low mass T Tauri stars, which contain co- 



tems which hav e flatter disks owin g to grain growth and 
settling. While Meeus et al. | ( 2001] ) invoked dust settling 
to explain the differences in flaring between Group I and II 
sources, it has also been suggested that the outer regions of 
some Herbig disks may experience self-shado wing, owing 
to a "puffed-up" inner circu mstellar disk rim (|Natta et al. | 
200T||Dullemond et al.|200T||Dullemond fc Dominik|2004a| 



pious amounts of circumstellar gas and dust ( Waters fc 
Waelkens|1998 1 . A wealth of observational evidence (IMan- 

■ tL if.- L i t nn>7 ll /">..j.„„::... t> I i nnn l I n ^ JL;„ „. 



nings & Sarg erit|1997| Oudmaijer fc Drew|1999| Mann ings 



fc Sargent||2000| |Vink et al.||2002| |Eisner et al.|2004| > sug- 
gests that the geometry of this gas and dust takes the form 
of a circumstellar disk. By contrast, our understanding of 
fundamental properties such as grain composition and size 
distribution, as well as the time evolution of these param- 
eters, are much less understood. 



Isella & Natta 2005 1. Recent models which include inflated 
inner disk rims have proven to be successful in explaining 
the near-IR a nd interferometric properties of several Her- 
big Ae stars (Ise lla et al.||2006"| ). Conclusive confirmation 
of the geometry of the outer disk regions of Herbig Ae/Be 
stars inferred by SED-based studies has generally not oc- 
curred. 

HD 163296 is a yo ung (4 Myr; 



van den Ancker et al 



2000[ ), nearby (122 pc; |van den Ancker et aL||1998[ ) Herbig" 
Ae star. While the star is not deeply embedded in a natal 



Analysis of ISO spectra led |Meeus et al. ( 2001[ ) to de- 
velop a 2-part evolutionary classification of Herbig sys- 
tems; Group I sources were characterized as being slightly 
younger systems with flared outer disks, while Group II 
sources were characterized as slightly more evolved sys- 



star formation dust cloud ( The et al.|1985 1, it still displays 
clear signs of active accretion via the presence of jets and 
Herbig Haro knots ( |Grady et al.||2000| [Devine et al.||2000 
Wassell et al. 120061). Hence, HD 163296 appears to be m a 



transition phase between optically thick, extremely young 
pre-main-sequence stars and the much more optically thin, 
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near zero-a 


ge-main-sequence debris disk type stars ( 


Ardila 


et al. 2004 


Kalas et al.|2005| 2006| Golimowski et al. 


2006). 



TT99M 

imaging ([Mannings fc Sargent||1997[ |G7ady et al.| [2000 
Isella et aT 2007p have confirmed the presence of a cir- 
cumsteiiar disk associated with HD 163296. HST STIS 
white-light coronagraphic observations spatially resolved 
the outer regions of this disk, and detected evidence of disk 
structure including an inner annulus of reduc ed scattering 
and a bright outer ring or ansae of material ( Grady et al. 
2000). The inner region of HD 163296's circumsteiiar disk, 
as diagnosed by near-IR (0.8 - 5 /im) SED monitoring, ex- 
hibits evidence of v ariability, possibly owing to changes in 
the inner disk wall ( |Sitko et al.||2008| ). 

Attempts to link processes which dominate the inner 
regions of Herbig Ae disks to morphological features ob- 
served in the outer disk regions, i.e. determining the rela- 
tive roles of self-shadowing, dus t settling, and disk flaring 
(Dullemond & Dominik 2004a I, requires a wealth of ob- 
servational data which diagnose the inner and outer disk 
regions. In this paper, we present a multi-epoch, multi- 
color view of the HD 163296 scattered light disk, provid- 
ing a significant improvement in the documented behavior 
of the outer regions of HD 163296's circumsteiiar disk. 
We describe the observational details and data reduction 
procedures applied to our ACS data in Section [2] In Sec- 
tion [3] we document the basic features of the HD 163296 
scattered light disk in our ACS data, and compare these 
results to earlier epoch scattered light imaging in Section 
[4] The possible mechanisms which might explain the ob- 
served scattered light variability are discussed in Section 
[5] We provide a summary of our main results in Section 

H 

2. OBSERVATIONS AND DATA REDUCTION 

We obtained coronagraphic imaging of HD 163296 in 
the HST ACS F435W filter via the cycle 12 GTO pro- 
gram 9987, and in the ACS F606W and F814W filters via 
the cycle 11 GTO program 9295. All observations utilized 
the coronagraphic mode of the High Resolution Channel 
(HRC) camera, along with its l'/8 diameter occulting spot, 
yielding a 29'.'2 x 26'.'2 field of view at a pixel scale of / /028 
x 07025. The FWHM of point spread function (PSF) of 
the HRC is 50 mas in V. Our science target data were ob- 
tained at two different HST roll angles to better facilitate 
the separation of PSF artifacts from real structures. The 
A3IV star HD 145570, selected for its spectral match to 
HD 163296, was used as a PSF reference, and was observed 
in orbits immediately following our HD 163296 exposures. 
A summary of the basic parameters of our dataset is pre- 
sented in Table [T] 

The standard HST pipeline was used to perform ba- 
sic initial data reduction procedures. We used the best 
available calibration files for this processing, including a 
correction for the different flat field illumination of the re- 
gion nearby the coronagraphic spot. CR-split images were 
combined using standard pipeline techniques. 

Techniques to identify and remove PSF signatures from 
HST ACS coronagraphic data are well documented in th e 
literature ( Clampin et al. 2003 Golimowski et al. 2006 1 . 
Following these established procedures, we first normalized 
and aligned the non-distortion corrected images of our PSF 



reference star to HD 163296 in an iterative manner, using a 
cubic convolution interpolation function. Residual align- 
ment errors in the re sultant best subtractio n are on the 
order of ±0.05 pixels ( |Pavlovsky et al.||2005[) and ou r nor- 
malization is accurate to ~2% (Clampin et al. 2003). Fol- 
lowing the subtraction of a fully registered and scaled PSF, 
we used MultiDrizzlc to correct our data for the known ge- 
ometric distortion of the HRC field of view and rotate all 
images into a common orientation. 

Next, SYNPHOT, the synthetic photometry package 
within STSDAS, was used to determine the correction fac- 
tor needed to calibrate our data to an absolute photomet- 
ric scale. As a template for HD 163296, we used a synthetic 
spectrum of 69 Her ( A2V) , normalized to the V-band mag- 
nitude of HD 163296, 6.899; a synthetic spectrum of 59 Her 
(A3III) normalized to a V-band magnitude of 4.940 served 
as a template for HD 145570. The ratio of the stellar flux 
between HD 163296 and HD 145570, extracted both from 
our coronagraphic images and our synthetic representa- 
tions of these data, were consistent to < 2% over multiple 
filter combinations, indicating our choice of synthetic spec- 
tra were reasonable. All PSF-subtracted, distortion cor- 
rected images were normalized to the synthetic flux of HD 
163296, and the different roll angle images were combined 
to increase the effective signal to noise ratios (SNR) . Note 
that all surface brightnesses reported in the native ACS 
filter bandpasses are based on the STMAG photometric 
system, and have been corrected for the known achromatic 
52.5% reduction of flux induced by the coronagraph. 



3. DISK CHARACTERISTICS DERIVED FROM HST ACS 
OBSERVATIONS 

Figures T] [2] and [3] present the PSF-subtracted, roll 
angle combined, calibrated images of HD 163296 in the 
F435W, F606W, and F814W filters respectively, using a 
square-root intensity stretch and a 3 pixel Gaussian ker- 
nel smoothing function. We have masked regions interior 
to a radial distance of 2'.'0, which are dominated by PSF- 
subtraction residuals. Note that several apparent features 
in these figures are actually well known instrumental ar- 
tifacts, whose ch aracteristics and origin are discussed in 
depth elsewhere dCla mpin et al. 2003 |Krist et al. 2005a 



Gonzaga et al.||2005p. All narrow radial spikes are residu- 
als due to incomplete PSF subtraction. Additionally, im- 
age imperfections arising from the 3'.'0 diameter occulting 
spot and its occulting finger are clearly seen in all images 
(WNW of the central source in the F435W filter, ESE of 
the central source in the F606W and F814W filters). This 
occulting finger and spot also projects "horn-shaped" ar- 
tifacts located at a distance of 2'.'25-3'.'0 from the central 
source. 

Three Herbig-Haro knots and a bipolar jet are clearly 
visible in our multi-epoch images; thes e features have been 
identified in previo us inv estigations dGrady et aT 2000 
Devine et aT||2000[ ), and |Wassell et al.| ( |2006p Has usee 
some of these data to estimate both the proper motion of 
the knots and the mass-loss rate of the HD 163296 system. 
Having identified these known features and artifacts of our 
data, we now focus on describing the radial and azimuthal 
surface brightness, color, and morphological features of the 
HD 163296 disk. 
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3.1. Disk Surface Brightness 

While the surface brightness profile of portions of the 
HD 163296 disk were c haracterized via wh ite-light STIS 
coronagraphic imaging (Grady et al. 20001, the large oc- 
culting wedges used by this coronagraph prevents an in- 
vestigation of the full radial and azimuthal structure of the 
disk. Our ACS images do not suffer from such artifacts, 
hence we are able to provide the first complete, multi-color 
assessment of the HD 163296 disk. 

We first calculate the disk surf ace brightnes s alon g the 
semi-major ax i s, as reported by Grady et al. (2000) and 



Wassell et al.| ( |2006| ). We extracted 10 pixel wi3e(0'.'25) 
radial profile cuts at position angles of 132° and 312° mea- 
sured north to east, and subsequentlycombined these data 
in 4 pixel (O'.'IO) bins (see Figure H}. As illustrated in 
Figure [5] for our F606W data, these radial profiles fol- 
lowed three distinct power-law trends. Between ~272 to 
~2'/8, evidence of scattered light is apparent; however, it 
is likely that PSF-subtraction residuals at least partially 
contribute to this measured flux. Clear evidence of the 
scattered light disk is seen in all filters from ^2'.'9 - ~4'.'4 
(350-540 AU) . The abrupt turnover in the profiles at radial 
distances >4'.'4 (540 AU) represents the termination of our 
detection of the disk; this boundary could cither indicate 
the physical outer radius of the disk's upper atmosphere 
or merely the radial location at which the disk becomes 
extremely optically thin. 

Table [2] lists the exponents of the best power law fits 
derived along the disk major axis in each of our three fil- 
ters. Concentrating on the ^2'.'9 - ~4'.'4 regions over which 
we are clearly seeing the scattered light disk, several fea- 
tures and trends are apparent. At the level of the SNR 
of our data, we do not detect clear evidence of differences 
between the surface brightness behavior in the SE versus 
NW quadrants of the disk in the F606W and F814W fil- 
ters. Note that the SE quadrant of our F435W disk major 
axis radial profile cut is "contaminated" by the presence 
of an ansa feature, which will be discussed in Section 3.3. 

3.1.1. Disk Size 

A large range of sizes h ave been suggested fo r the HD 
163296 circumstellar disk. Dominik et al. (20031 assumed 
an o uter radius of 50 A U m their SED fitting of the sys- 
tem. Isella et al. (2007) reported that CO emission, likely 
produced in the upper layers of the disk, extended out 
to a radial distance of 540 AU; moreover, they reported 
an unexpected deficit of mm continuum emission beyond 
200 AU which they interpreted as evidence that the disk 
was s t rongly depleted of dust outside of 200 AU. |Isella| 
et al. (20071 speculated that this deficit could be due to 
the inward migration of large bodies in the disk or due to 
clearing by a low mass companion. 

Our spatially resolved imagery reveals that light scat- 
tering off of dust residing in the upper layers of the HD 
163296 disk extends out to a radial distance of ~540 AU. 
Given our detection, we suggest that there is no physical 
gap of dust at radial distances g reater than 200 AU as hy- 
pothesized by llsella et al. (20071. Rather, we suggest that 



so me of t he alternate interpretations discussed in |Isella et] 
al. ( 2007 1 might provide more feasible explanations for the 
non-detection of dust mm continuum emission beyond 200 
AU, including grain growth and/or a change in the disk 



geometry beyond 200 AU (assumed by Isella et al.|2007 to 
be a fully flared disk). 

3.1.2. Disk Geometry 

Power-law fits to scattered light imagery provide one 
means to diagnos e the geometry of H erbig Ae circumstel- 



lar disks (see e.g. Grady et al 



20071. As seen in Table 2 



the surface brightness of the HD 163296 disk in the ACS 
F435W, F606W, and F814W filters follow a ~r~ 4 depen- 
dence. These data clearly deviate from the expected be- 
havior of a continuously flared disk, which would follow a 
' 2 dependence, or a geom etrically thin disk, which would 



follow a r 3 depend ence (Whitney & Hartmann 19921. 
iDominik et al.l ( 2003 ) also reached the conclusion that the 



HD 163296 disk was not continuously flared based on SED 
modeling, though some of the parameters adopted in their 
model which produced this result (i.e. a disk outer radius 
of 50 AU) are dubious. 

3.2. Disk Color 

Our multi-filter HST ACS dataset also allows us to pro- 
vide the first measurement of the color of HD 163296's 
scattered light disk. We first calculated the surface bright- 
ness at radial distances of 3'.'0 and 3'.'5 along the disk semi- 
major axis (PAs of 132° and 312°). The median flux within 
l'.'O x l'.'O apertures were determined at each location, and 
complimentary measurements (PAs of 132 and 312) were 
averaged to produce the surface brightnesses listed in Ta- 
ble [3] Given our previous assessment that PSF-subtraction 
residuals could be contributing to the observed scattered 
light inside a radial distance of 2 / .'9, we caution that our 
aperture photometry at 3'.'0 similarly might be marginally 
contaminated: by contrast, the aperture photometry at 
3'.' 5 should be free from such effects. Background sky sur- 
face brightnesses (Table |3| were determined by averaging 
the median pixel values in two l'.'O x l'.'O apertures located 
9'.'0 from the central star. The standard deviation of these 
background sky apertures were used to provide an esti- 
mate of the photometric uncertainties within our science 
apertures at 3'.'0 and 3'.' 5 along the disk semi-major axis. 

The surface brightnesses in each of our ACS filters (Ta- 
ble [3j) were converted to the standard BVI system by iter- 
ation of the transformation equation outlined by |Sirianni] 
et ai.|fl2005| ); TMag = SMag + cO + cl * TC + c2 * TC Z , 
where TMag is the Johnson B, V, or I magnitude, SMag 
is the ACS filter magnitude in the OBMag photometric 
system, TC is the OBMag color difference of the filters 
used, and cO, cl, and c2 are transformation coefficients 



tabulated by Sirianni et al. (2005|. Note that the B-band 
surface brightness and (B-V) color were calculated using 
multi-epoch data, while V- and I-band surface brightnesses 
and (V-I) color were calculated using single epoch data 
(see Table [I]). 

The observed disk (V-I) color 3'.'5 from the central star 
along the semi major axis, 1.1, is significantly redder than 
literature measure ments of the stellar (V-I) color, 0.15 
To better assess the global behav- 



(The et al. 



1985). 



ior of the disk's color, we transformed our F606W and 
F814W imagery to the Johnson V- and I-band filters us- 
ing the aforementioned techniques. The resultant (V-I) 
image of the scattered light disk, plotted on a linear scale 
and smoothed with an 8 pixel Gaussian kernel, is shown in 
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Figure [6] From inspection of Figure [6] we conclude that 
our data provide no clear evidence of a spatial dependence 

in the (V-I) color of the scattered light disk. 

With the exception of AU Mic (|Krist et al.||2005a[ ), the 
)lors of many other resolved disk systems are similarly 



colors of many other resolved disk systems arc similarly 
redde r than their h ost stars (HR 4796A , Schneider et al. 
1999j HP 141 569A, |Clampin et aT1[2003l GG T auri, |lgist 



et al. 2005b f3 Pic, Golimowski et al. 20061. As dis 



cussed for these aforementioned systems, disk colors have 
been used to provide tentative constraints on the relative 
size distribution and/or chemistry of grain populations, al- 
though degeneracies clearly exist in such an analysis (e.g. 
McCabe et al.|2002[ ) . One possible interpretation of a red 
disk color is that some amount of grain processing has 
occurred (i.e. grain growth; chemistry evolution). Grain 
growth and grain settling is a ph enomenon which is be- 
lieved to occur in Herbig Ae disks (|van Boekel et alT|2 003 
Acke et alT||2004| |Pullemond fc Pominik||2004a|b| , and is 
generally associated with systems whose outer disks are 
not continuously flared. Previous analysis of the IR SEP 
of HP 163296 has led to suggestions that thi s disk system 
has e xperienced grain growth and settling (Meeus et al. 
2001 1 . The detection of a non-continuously flared disk and 
red (V-I) disk color from our spatially resolved scattered 
light imagery is thus consistent (though not conclusive) 
with such phenomenon occurring. We also note that our 
HP 163296 imagery does not exhibit significant evidence 
of a scattering asymmetry due to differences in forward 
and back scattering of grai ns in the disk, as i s seen in 
systems such as GG Tauri (Krist et al. 2005b I. Future 



modeli ng efforts which uti lize both the known SEP prop 



erties (Meeus et al. 20011 and scattered light properties 
(the present study) of the HP 163296 disk would help to 
better constrain the grain size distribution in a more quan- 
titative manner. 

3.3. Morphological Features 

The superlative spatial resolution of the HST ACS al- 
lows for the search of both radial and azimuthal morpho- 
logical structures within the HP 163296 scattered light 
disk. Simple visual inspection of our F606W and F814W 
images (Figure [7]) reveals no clear evidence of disk struc- 
ture. By contrast, our later epoch F435W data (Figure [71 
clearly exhibits evidence of an ansa-like structure in the 
southeast (SE) quadrant of the disk, located at the same 
position of one of two ansae noted in archival HST STIS 
coronagraphic observations ( Grady et al.|2000 l. The ansa 
is visible and stationary in both of the roll angle images 
at F435W, whereas PSF-subtraction artifacts rotate with 
the roll angle; moreover, this structure is not associated 
with the "horn-shaped" artifacts attributable to the ACS 
large occulting spot and finger. 

We extracted 10 pixel wide (CK'25) median averaged ra- 
dial profile cuts from our data to parameterize this disk 
structure. The bottom panels of Figure [8] depict the ra- 
dial surface brightness profiles, in 3 pixel (0'.'075) radial 
bins, at position angles of 140° N to E (SE disk) and 320° 
N to E (NW disk) for our three band-passes. At its bright- 
est point, a radial distance of ~3'.'l-3'.'2, the structure has 
a F435W surface brightness of — 19.5 mag arcsec -2 , and 
is preceded by sudden dimming in surface brightness. At 
slightly different position angles, 130° N to E and 310° 



N to E (see top panel, Figure [8]) , the feature is slightly 
brighter (F435W —19.0 mag arcsec -2 ), characterized by 
a less steep dimming in surface brightness, and reaches a 
peak flux at a slightly more interior radial distance. At 
both position angles, no corresponding structures are seen 
cither in the NW cut of the F435W data or the SE and 
NW cuts of the F606W and F814W data. Note that the 
apparent "feature" in the F814W filter (and less apparent 
in the F606W filter) data at 3'.'5 in the NW disk quad- 
rant is merely a background star. It is quite clear that the 
overall surface brightness of the F435W data at the radial 
distance of the feature (— 3'.'0) is brighter in the SE region 
of the disk than in NW region of the disk. 

4. THE MULTI-EPOCH BEHAVIOR OF THE HD 163296 
SCATTERED LIGHT DISK 



Grady et al. 



^ ([2000) presented white-light HST STIS 
coronagraphic observations of HP 163296, obtained in 
September 1998. The STIS observations were character- 
ized by a featureless signal from l'/5 to 2'.'5, a darker zone 
at 2765 of width Off 4, and two bright ansae peaking at 
2'.' 9 in the SE quadrant and 3 " 2 in the NW quadrant (see 
Figure [7]). Grady et al. (20001 considered two primary ex- 
planations fOTtn^^b^erve^~Bright ansae and dark lanes in 
their STIS imagery, changes in the dust grain composition 
and dynamical clearing by planetary-like bodies, although 
they noted other possibilities could not be excluded. 

Comparison of these archival white-light (X ce ntrai — 
5850 A, AA = 4410 A) STIS data with our multi-epoch 
F435W (X cen t r ai = 4297 A, AA = 1038 A), F606W 
(Xcentrai = 5907 A, A A = 2342 A), and F814W (A celltrai 
= 8333 A, AA = 2511 A) ACS observations requires cau- 
tion, given the different band-passes sampled by each indi- 
vidual observation. Our single-epoch F606W and F814W 
observations span most of the STIS bandpass. To serve 
as a crude, approximate template for comparison to the 
archival STIS images, we have combined our ACS F606W 
and F814W data after first normalizing the images to ac- 
count for the different throughput of these filters, —0.16 for 
F814W and -0.26 for F606W. Hereafter we refer to this 
combination as our ACS F606W+F814W dataset. Simi- 
larly, under the assumption that the HD 163296 disk sur- 
face brightness is not variable, we have combined our dif- 
ferent epochs of ACS data, again after normalizing the 
images to account for the different throughput of the fil- 
ters, -0.16 for F814W, -0.26 for F606W, and -0.19 for 
F435W. We will refer to this latter combination, which 
covers the majority of the STIS clear filter bandpass, as 
our ACS F435W+F606W+F814W dataset. 

4.1. Disk Surface Brightness 



Mirroring the analysis performed in Section 3.2 we cal- 
culated the average disk surface brightness along the semi- 
major axis of the disk, at radial distances of 3'.'0 and 3'.'5 
from the central star. These surface brightness measure- 
ments are listed in Table [3J we also report the surface 
brightness of the STIS epoch data in the Johnson V-band 
in Table [3] We stress that this quoted V-band surface 
brightness, calculated via a transformation equation which 
incorporated our observed (V-I) disk color, is only a crude 
approximation. 
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The 1998 STIS epoch surface brightness along the major 
axis of the HD 163296 scattered light disk is significantly 
fainter than all other measurements in all other filter band- 
passes. In the Johnson V-band system, the 1998 STIS ob- 
servation at a radial distance of 3'.'5 from the central star 
(20.1 mag arcsec -2 ) is ~1 mag arcsec -2 fainter than the 
2003 ACS V-band surface brightness (19.2 mag arcsec -2 ; 
see Table [4}. In the instrumental filter system, the 1998 
STIS clear filter surface brightness at a radial distance of 
3'.'5 from the central star (20.3 mag arcsec -2 ) is similarly 
~1.0 mag arcsec -2 fainter than that measured in both of 
our ACS "white-light" filter approximation datasets (19.3 
and 19.1 mag arcsec -2 respectively; see Table [3]). As the 
1998 epoch STIS data are significantly fainter than that 
measured in each bandpass of our ACS data, it is clear 
that no combination of the ACS data can reproduce the 
behavior of the STIS scattered light surface brightnesses. 
We suggest that the observed multi-epoch behavior of the 
disk's surface brightness is one piece of evidence which 
suggests that the system's scattered light disk is variable. 

4.2. Morphological Features 



H 



Grady et al. ( 2000 ) reported the presence of two ansae in 
the white-light STIS coronagraphic image of HD 163296 
obtained in 1998, with the brightest ansa located in the 
NW quadrant of the disk. In Section |3.3| we reported a 
null detection of any morphological features in our 2003 
epoch data (F606W and F814W filters) and the detection 
of a single ansa structure in the SE quadrant of our HST 
ACS F435W data. The variability in the visibility of these 
morphological features arc easily seen in Figure [7] 

To compare and contrast these features in a quantita- 
tive manner, we have followed the technique outlined in 
Section 3.3 and extracted radial surface brightness profiles 
for the ACS F606W+F814W and STIS data (see Figure 

t. Note that our ACS radial profile cuts correspond to 
e median average of a 10 pixel wide (0 / .'25 arcsec) aper- 
ture which were subsequently averaged in 3 (radial) pixel 
bins, while our STIS radial profile cuts correspond to the 
median average of a 5 pixel wide (0f.'25 arcsec) aperture 
which were subsequently averaged in 3 (radial) pixel bins. 

It is clear that the ansa structure we reported in our 
ACS F435W data in the SE quadrant of the disk (top left 
panel, Figure |9| occurs at the same r adial distance as th e 
ansa in the SE quadrant reported by Grady et al. (20001; 
hence it is plausible to suggest that these morphological 
features are in fact the same phenomenon. Interestingly, 
while the ansa is brighter in the NW quadrant than in 
the SE quadrant in the STIS data, a nearly opposite sce- 
nario characterizes the ACS F435W data, such that the 
ansa is present in the SE quadrant and absent in the NW 
quadrant. The ansa features which characterize the HD 
163296 scattered light disk not only exhibit variability in 
the epochs in which they are discernible, but also exhibit 
different surface brightnesses in the epochs in which they 
are visible. 

The non-detection of these morphological features in our 
2003 epoch data are unlikely to be due to the sensitiv- 
ity limits of our data. Ignoring the obvious differences in 
band-passes, the observed surface brightness of our ACS 
data which fail to detect any morphological features (2003; 
F606W, F814W) are actually brighter than that observed 



in our ACS data in which structure is detected (2004; 
F435W). Moreover, the limiting magnitude of our 2003 
ACS data, depicted as dashed horizontal lines in Figure 
is one magnitude arcsec -2 below the measured surface 
brightness of the morphological feature in our 2004 epoch 
F435W data, inside of a radial distance of 3'.'5. 

One might also reason that, given the different band- 
passes probed by our multi-epoch ACS and STIS data, 
perhaps the visibility of the morphological structure is an 
optical depth effect, and is extremely color dependent (i.e. 
produced by small particle size scatterers) . Given the very 
red (V-I) color of the disk (Section 3.2 1, we suggest that 
this scenario is dubious. Moreover, it is difficult to envision 
how an optical depth phenomenon could simultaneously 
explain both the visibility pattern (visible in 1998 white- 
light; 2004 F435W) and surface brightness variability (NW 
ansa brighter than SE ansa in 1998; SE ansa brighter than 
(non-detected) NW ansa in 2004) of these morphological 
features. Rather, we suggest that the visibility and surface 
brightness variability of these morphological features are 
two additional pieces of evidence which suggest that the 
system's scattered light disk is variable. 

We also note that unresolved linear spectropolarimet- 
ric monitoring of HD 163296 has revealed evidence of a 
time variable intrinsic polarization component (Bjorkman 
2007, personal communication). Such data indicate vari- 
ability in the net (unresolved) scattered light behavior of 
the system. As the integrated scattered light of HD 163296 
exhibits variability, we suggest it is perhaps not surprising 
to find evidence that when this scattered light is resolved 
via coronagraphic imaging, that one finds multiple avenues 
of evidence indicating variability. 

5. ORIGIN OF THE OBSERVED SCATTERED LIGHT 
VARIABILITY 

Our analysis of multi-epoch, multi-color spatially re- 
solved imagery of HD 163296 suggests that its scattered 
light disk is variable; we observe a ~1 mag arcsec -2 change 
in the overall surface brightness of the disk between epochs 
of observations, variability in the visibility of ansa struc- 
ture (s), and azimuthally asymmetric variability of the sur- 
face brightness of these ansa(e) , when they are visible. We 
now consider the possible origins of this variability, in the 
context of what is already known about the nature of the 
HD 163296 disk. 

Previous interpretations of the bright ansae and dark 
zones which chara cterize the STIS scattered light data 
( |Grady et aL|2000| outlined numerous mechanisms which 
could explain the origin of the observed phenomenon, in- 
cluding dynamical clearing, grain composition and/or size 
distribution differences, changes in the degree of disk flar- 
in g , and c hanges in the number density of grains. | Grady et| 
al. (2000) also postulated that, if the darker zones were in 
fact a partially cleared annulus, then a ~0.4 Mj Mp j ter gas 
giant planet might be responsible for producing the disk 
clearing. While our 2003 epoch ACS F606W and F814W 
data clearly detected the HD 163296 scattered light disk, 
no evidence of any dark lanes were observed. Thus it ap- 
pears unlikely that the m orphological "dark lane" identi- 



fied by Grady et al. ( 2000 1 corresponds to a cleared region 
of the disk; we suggest that the scattered light data present 
no clear evidence of the presence of a planetary mass body. 
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We consider two classes of mechanisms which could 
induce variability in our scattered light observations, 
changes in the fundamental properties of the scatterers 
and changes in the illumination of these scatterers. Our 
data indicate that the scattered light disk is variable on at 
least a time-scale of several years, which strongly suggests 
that we can exclude changes in the fundamental properties 
of the scatterers (i.e. grain composition, grain size distri- 
bution, grain number density, and degree of disk flaring) 
as a plausible mechanism. 

Changes in the illumination of the disk, from the per- 
spective of the scatterers, could include either (or both) 
stellar variability or variations in the number of stellar 
photons which reach the scatterers. We are not aware of 
any evidence which suggests that HD 16329 6 experiences 
UX-Ori-like photometric ou tbursts (Waters & Waelkens 
1998 Dullemond et al. 2003); rather, the star appears to 
have a generally stable photometric brightness of V=6.88, 
with a maxi mum variability range (high to low) of 0.42 
magnitudes ( Herbst fc Shevchenko||1999 l. Direct imaging 
of HD 163296 was obtained in each of our ACS filters at the 
epoch of our coronagraphic observations; however, these 
data were heavily saturated and did not use a gain set- 
ting which sampled the full well depth of the HRC, which 
would have permitted accurate phot ometry to be extracted 
( |Gilliland|2004l|Sirianni et al.|2005| . Thus, we do not have 
any observational evidence which could link HD 163296's 
observed scattered light variability with stellar variabil- 
ity. Moreover, it is not clear that brightening/dimming 
of the central star could explain the visibility and surface 
brightness variability patterns observed in the disk ansae 
features. Rather, we suggest that HD 163296's scattered 
light variability is related to changes in the number of pho- 
tons which reach the disk scatterers. 

Self-shadowing of the outer regions of Herbig Ae disks 
is a phenomenon predicted by theory (Natta et al.|2001 



Dull emond et al.|2001| |Dullemond fc Dominik|2004a||TserIa 



fc Natta||2Q05 ) , and variable illumination of e nvelopes has 
been r eported in observations of GG Tau flKrist et al~ 
2005b| ) and Hubble's Variable Nebula QLightfoot ||1989|) 
One suggested method to produce a self-shadowed disk 
is by inflation of the inner disk wall, which prevents co- 
pious amounts of phot ons from reaching at least part of 
the outer (flatter) disk. Sitko et al. (20081 has noted that 
the infrared spectral energy distribution (IR SED) of HD 
163296 can vary by ~30% on time-scales of 3 years or 
less, which they interpret as evidence of structural changes 
in the region of the disk near the dust sublimation zone. 
The mechanism driving this inner disk variability is not 
known, although Sitko et al. (2008) discuss several possi- 
bilities: a-disk instabilities, magneto-rotational instabili- 
ties, an X-wind origin, and perturbations by a planetary- 
mass body. Given the non-continuously flared geometry 
of the HD 163296 outer disk, as inferred from our ACS 



scattered light data (Section 3.1 1, it is plausible that even 
small-scale changes in the scale height of the inner disk wall 
could induce a variable amount of self-shadowing of the 
outer disk. Thus we suggest that variable self-shadowing is 
one plausible mechanism which might explain the observed 
variability of the outer scattered light disk. Coordinated, 
contemporaneous observations of the inner and outer disk 
regions of HD 163296, namely follow-up IR SED and coro- 



nagraphic monitoring, would provide a robust means of 
investigating this proposed mechanism. 

Recall that the bright ansa(e) features which sometimes 
characterize the HD 163296 scattered light disk were only 
observed during epochs in which the overall disk surface 
brightness was low (1998, 2004; see Figure [9] and Tables 
[3]and|4|. Thus, in our working picture of the HD 163296 
scattered light disk, with admittedly few data points to 
work with, these morphological features are only detected 
during periods of enhanced self-shadowing. We suggest 
that the source of these observed bright features is not a 
density enhancement of scatterers (i.e. a clump), as one 
would expect such a clump would contribute a localized 
enhancement of scattered light not only during enhanced 
self-shadowing periods, but also during epochs character- 
ized by less self-shadowing. Rather, we speculate that 
these bright features signify a change in the spatial distri- 
bution of the scatterers, i.e. the ansa(e) represent a region 
in the disk which is modestly less flattened than other 
surrounding locations. During periods of increased self- 
shadowing these structures would partially escaped the 
shadow created by the inner disk wall. During periods 
of less self-shadowing, this region would be illuminated 
similar to its neighboring disk regions, thus as long as 
the structures do not have a significantly higher density 
of scatterers, the morphological feature would not be con- 
spicuous in scattered light profiles. While the origin of 
the perturbations which give rise to these localized fea- 
tures is unknown, structure in the tenuous upper layers 
of protoplane tary disks is predicted to o ccur by theory. 
For example, Jang-Condell & Boss (20071 predicted such 
a phenomenon would occur in the inner regions of proto- 
planetary disks as a byproduct of planet formation via the 
disk instability model. 

Finally, we consider the origin of the surface brightness 
variability of the ansa(e) struct ures, during the epochs in 
which they were visible (Section 4.2 Figure [9]). If the vis- 
ibility of these features is linked to the behavior of the 
inner disk wall, as we have suggested, then it is likely that 
the surface brightness variability of these features is also 
linked to events in the inner disk wall region. We speculate 
that the inner disk wall might be azimuthally asymmet- 
ric, hence cast an azimuthally asymmetric shadow on the 
outer disk. Given the orbital period of material located 
at the inner wall sublimation radius, <1 week (Sitko 2007, 
personal communication) , the time-scale for such a asym- 
metry to produce a brighter NW versus SE illumination 
(and vice- versa) is much less than the maximum time-scale 
of the observed ansa variability, several years. We suggest 
that a series of scattered light images of HD163296, ob- 
tained in time steps of ~1 to several weeks and during a 
period of enhanced self-shadowing, would provide a direct 
test of our interpretation of the system. 

6. SUMMARY 

We have presented the first multi-color, multi-epoch 
analysis of resolved optical scattered light imaging of the 
Herbig Ae star HD 163296. To summarize the observa- 
tional properties of these data: 

1. We spatially resolved the HD 163296 scattered light 
disk over radial distances of 2 / .'9-4'.'4 (350-540 AU) 
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in the HST/ACS F435W, F606W, and F814W fil- 
ters. 

2. Radial profiles of the surface brightness along the 
semi-major axis of the disk follow a ~r~ 4 power law 
behavior, which is indicative of a non-continuously 
flared disk. 

3. The (V-I) color of the disk at 3'.'5, ~1.1, is signifi- 
cantly redder than the stellar (V-I) color, 0.15. This 
red disk color appears to be spatially uniform at 
the SNR of our data, and is consistent with that 
expected for a disk which has experienced grain 
growth and is at least partially self-shadowed. 

4. A single ansa structure is present in the SE disk 
quadrant of the 2004 epoch ACS data (F435W fil- 
ter), while no structure was observed in the 2003 
epoch ACS data (F606W and F814W filters), de- 
spite the fact that the limiting detection magnitude 
of these latter data was sufficient to probe structure 
at the surface brightness observed in 2004. This 
morphological feature is spatially coincident with 
the fainter ansa rep orted in 1998 STIS observations 
( |Grady et al.|2000[ ) of HD 163296; the brighter ansa 
seen in the NW disk quadrant of the STIS data are 
absent from our ACS epoch data. 

5. The scattered light disk was observed to be signifi- 
cantly brighter in each of three filters of ACS obser- 
vations from 2003-2004, as compared to white-light 
STIS observations in 1998. Combining the multi- 
color ACS data to crudely approximate the STIS 
bandpass suggests that the ACS epoch data are ^1 
mag arcsec -2 brighter than the STIS data. Along 
with the variability of the visibility and surface 
brightness of ansa structure(s) in the disk, these 
results suggest the scattered light disk of HD 
163296 is variable. 

We expect that the characterization of the basic behav- 
ior of HD 163296's scattered light disk that we have pro- 
vided in this study will serve as an important foundation 
to future efforts to model the multi-wavelength behavior 
of the system. Although speculative, we suggest that one 



plausible explanation for the origin of the observed vari- 
ability is: 

1. The scale height of the inner disk wall is believed 
to inflate and deflate on time-scales of less than a 



fe w year s, based on IR SED monitoring (|Sitko et| 
al. 20081. The variable scale height of the inner 
disk wail could induce variable self-shadowing of the 
outer disk, hence produce the observed overall vari- 
ability of the scattered light disk. 

2. We suggest that the ansa structure (s), which ap- 
pear to be discernible only during periods of en- 
hanced self-shadowing, represent a localized region 
of scatterers which are at a scale height (at least 
marginally) above the projected shadow. As the 
features are not discernible during periods of more 
complete illumination of the outer disk, we believe 
that a localized enhancement in the scale height of 
the disk is more likely to produce the observed phe- 
nomenon than a localized density enhancement of 
the scatterers (i.e. a clump). 

3. The relative surface brightness variability of the 
NW versus SE ansa structures during epochs in 
which they are visible suggests that the inner disk 
wall might be azimuthally asymmetric, hence pro- 
duce an azimuthally asymmetric shadow on the 
outer disk. Based on the short (<1 week, Sitko 
2007 personal communication) orbital time-scale of 
material located at the inner wall, we suggest that 
a series of resolved scattered light images of HD 
163296, obtained during an epoch of enhanced self- 
shadowing, would provide a test of this particular- 
suggested phenomenon. 
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Table 1 
Summary of Observations 



Object 


Date 


Instrument 


Filter 


Expos. Time 


PA 


Comment 


HD 


163296 


16 Aug 2004 


HST ACS HRC 


435W 


2275 sec 


281 


roll angle 1 


HD 


163296 


16 Aug 2004 


HST ACS HRC 


435W 


2480 sec 


251 


roll angle 2 


HD 


145570 


16 Aug 2004 


HST ACS HRC 


435W 


2320 sec 


280 


PSF star 


HD 


163296 


25 Mar 2003 


HST ACS HRC 


606W 


2350 sec 


76 


roll angle 1 


HD 


163296 


25 Mar 2003 


HST ACS HRC 


606W 


2350 sec 


104 


roll angle 2 


HD 


145570 


25 Mar 2003 


HST ACS HRC 


606W 


900 sec 


94 


PSF star 


HD 


163296 


25 Mar 2003 


HST ACS HRC 


814W 


2100 sec 


76 


roll angle 1 


HD 


163296 


25 Mar 2003 


HST ACS HRC 


814W 


2260 sec 


104 


roll angle 2 


HD 


145570 


25 Mar 2003 


HST ACS HRC 


814W 


900 sec 


94 


PSF star 



Note. - 



Note that the PA listed in column six refers to the position angle of the V3-axis (roll-angle) of the HST. 



10 



Table 2 

Power Law Exponents of Disk Major Axis Radial Profiles 



Filter-Disk Quadrant Region 1 Region 2 Region 3 
2'.'2-2'.'8 2'.'9-4'.'4 4'.'5-6'.'l 



F435W-NW 


1.4 


±1.1 


3.7 


±0.3 


2.2 


±0.6 


F435W-SE 


4.9 


±1.1 


5.3 


±0.4 


1.5 


± 0.5 


F606W-NW 


0.5 


±0.8 


4.6 


±0.2 


1.1 


±0.2 


F606W-SE 


-1.0 


±0.4 


3.9 


±0.1 


1.3 


±0.2 


F814W-NW 


2.4 


±0.5 


3.7 


±0.3 


0.4 


±0.3 


F814W-SE 


-0.7 


±0.6 


3.2 1 


± 0.4 


0.7 2 


±0.3 



Note. — Power law fits to the surface brightness along the south-east (Left) and north-west (Right) regions of the semi-major axis in 
the F606W filter. The radial profiles were determined from the median of 10 pixel-wide (0'.'25) cuts across the disk major axis, which were 
subsequently combined in 4 pixel (0"10) radial bins. Radial profile cuts at a disk position angle of 132°, measured north to east, sampled the 
southeast (SE) disk quadrant while cuts at a disk position angle of 312° sampled the northwest (NW) disk quadrant. Fit regions which were 
truncated from that listed in columns 2-4 include: 1 3"1 - 3"7, 4" 4 - 4"5; 2 4'.'6 - 6'.'1. 
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Table 3 
Surface Brightnesses 



Filter 


Disk Major Axis 
3'.'0 

(mag arcsec - 2 ) 


Disk Major Axis 
3"5 

(mag arcsec - 2 ) 


Background Sky S.B. 
9'.'0 

(mag arcsec" 2 ) 


ACS F435W 


19.6±0.1 


20.1±0.2 


22.3 


ACS F606W 


18.6±0.2 


19.2±0.2 


20.6 


ACS F814W 


19.1±0.2 


19.4±0.2 


20.8 


ACS (F606+F814) a „ e 


18.9±0.2 


19.3±0.2 


20.7 


ACS (F435+F606+F814) a „ e 


18.7±0.2 


19.1±0.2 


20.7 


STIS Clear 


20.0±0.1 


20.3±0.2 


21.8 



Note. — Disk surface brightnesses, calculated by median averaging l'.'O x l'.'O apertures placed at radial distances of 3'.'0 and 3"5 along 
the disk major axis (PA = 132°, 312° measured north to east) are cited in the STMAG photometric reference frame. The background sky 
surface brightnesses were determined by averaging the median pixel values within two l'.'O x 1"0 apertures located 9'.'0 from the central star. As 
discussed in Section]!] the ACS (F606W+F814W) at , e and ACS (F435W+F606W+F814W) o „ e data correspond to weighted averages of single- 
and multi-epoch ACs data respectively which were combined to crudely approximate the STIS clear filter bandpass. 
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Table 4 

Johnson Filter Disk Major Axis Surface Brightnesses and Colors 



Filter 


at/o 


3" 5 




(mag arcsec -2 ) 


(mag arcsec -2 ) 


B (ACS 2004) 1 


20.2±0.1 


20.7±0.2 


V (ACS 2003) 


18.5±0.2 


19.2±0.2 


I (ACS 2003) 


17.8±0.2 


18.1±0.2 


V (STIS 1998) 


19.6±0.1 


20.1±0.2 


(V-I) 


0.7±0.3 


1.1±0.3 



Note. — The disk surface brightnesses quoted in Table [3] were converted to the sta ndard Johnson UBVRI system. Transformations of the 
ACS data were achieved using the iterative technique outlined by |Sirianni et al. |2005 I . Note that all data characterized by a 1 reference were 
derived via multi-epoch ACS observations. 
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0.0 0.12 

Fig. 1.— The PSF subtracted, 2 roll angle combined image of HD 163296 in the F435W filter. The image, which has a field of view of 20'.'0 
x 20'.'0, has been plotted using a square root scale, smoothed with a 3 pixel Gaussian kernel. Regions inside a radial distance of 2'.'0 were 
dominated by PSF-subtraction residuals and were masked from the displayed data. 
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15 




Fig. 3. — HD 163296 in the F814W filter, using the same parameters as in Figurc[T] 



1G 




23 i .... i ... I 

-5 -4 -3 3 4 5 

Radius (arcsec) 

Fig. 4. — Ten pixel wide (0'.'25) median radial surface brightness cuts along the semi-major axis (PA = 132° and 312°) of the HD 163296 
scattered light disk. The data were subsequently combined into 4 pixel (O'.'IO) wide bins. The colored horizontal lines denote the background 
sky surface brightness in each filter. 
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-0.4 



-0.6 



-0.8 



-1.2 



-1.4 
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Log Radius (arcsec) 



Fig. 5. — Power law fits to the surface brightness along the south-east (Left) and north-west (Right) regions of the semi-major axis in 
the F606W filter. The radial profiles were determined from the median of 10 pixel-wide (0'/25) cuts across the disk major axis, which were 
subsequently combined in 4 pixel (0"10) radial bins. Three distinct fitting regions, 2'.'2-2'.'8, 2'/9-4'.'4, and 4'.'5-6'.'l, denoted in this figure by 
different color schemes, yielded the power law parameters listed in Table [2] A similar analysis was performed for the F435W and F814W 
data. 
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0.6 1.8 



Fig. 6. — A 14'.'5 x 14'.'5 (V-I) color image of the HD 163296 scattered light disk plotted on a linear scale, smoothed with a 8 pixel Gaussian 
kernel. Regions inside a radial distance of 2'.'0 are masked, and the large white circle corresponds to a radial distance of 4'.'0. The color of the 
scattered light disk is similar to that of the 3 H-H knots visible, and exhibits no significant evidence of systematic spatial variations at the 
SNR level of our data. 
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Fig. 7.— 1998 epoch HST STIS clear filter (upper left), 2003 epoch HST ACS F606W (upper right), 2003 epoch HST ACS F814W (lower 
left), and 2004 epoch HST ACS F435W (lower right) coronagraphic observations of HD 163296. The field of view for all images is 12'.' x 
12"0; all images are plotted on a square root stretch scale smoothed by a second order Gaussian function, except for the STIS image which is 
smoothed with a first order Gaussian. Regions inside a radial distance of 2'.'0 are masked, and the large black circle corresponds to a radial 
distance of 4''0. 
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Fig. 8. — Bottom Panels The median average of ten pixel wide (0'.'25) radial profile cuts, subsequently averaged in bins 3 radial pixel 
wide, were calculated at position angles of 140° and 320°, measured N to E. These cuts clearly depict the behavior of the spiral arm structure 
observed in the SE disk quadrant of the F435W data (at a radial distance of ~3'.'0-3''5). The dashed horizontal lines represent the limiting 
surface brightness magnitudes for each bandpass, while the vertical lines depict the three regions characterized by distinctive power laws in 
Table m Top Panels. The median average of ten pixel wide radial profile cuts, subsequently average in 3 (radial) pixel bins, calculated at 
position angles of 130°, 140°, 310°, and 320° for the F435W data. 
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Fig. 9. — Bottom Panels The median average of 0'.'25 arcsec wide radial profile cuts are shown for all datasets, calculated at position angles 
of 140° and 320°, and subsequently averaged in bins 3 radial pixels wide. Note the epochs characterized by low overall surface brightness 
profiles (STIS - pink, ACS F435W - blue) are also the epochs in which spiral arm structure(s) are observed in the disk. The dashed horizontal 
lines represent the limiting surface brightness magnitudes for each bandpass, while the vertical lines depict the three regions characterized by 
distinctive power laws in Table [2] Top Panels. The median average of 0''25 arcsec wide radial profile cuts are shown for all epochs which 
exhibit spiral arm structure(s). 



